home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / modellib / default_init.c < prev    next >
Encoding:
C/C++ Source or Header  |  1990-01-25  |  2.4 KB  |  82 lines
  1. /*
  2. ### Initialize all default parameters of a dynamical system to be installed ###
  3. */
  4.  
  5. int default_init()
  6. {
  7.     int i;
  8.     extern int ir,var_cur_index[2],param_cur_index[2];
  9.     extern int var_dim,param_dim,full_dim,func_dim;
  10.     extern int mapping_on,inverse_on,fderiv_on,enable_polar,enable_period;
  11.     extern int *aux_on,*aux_win_mode,*aux_rescale,**aux_cur_index,aux_max;
  12.     extern double *param,*param_min,*param_max,*func,*func_min,*func_max;
  13.     extern double *var_i,*var_polar_i,*var_min,*var_max,*var_polar_min,*var_polar_max;
  14.     extern double *period_len,pi,*xoffset,*int_yscal;
  15.     extern char *title_label,**var_label,**var_polar_label,**param_label,**func_label;
  16.  
  17.     /*periodicity of phase space variables (DIM=var_dim)*/
  18.     for(i=0;i<full_dim+param_dim;i++) period_len[i] =0;
  19.     /* primary phase space variable label (DIM=var_dim)*/
  20.     for(i=0;i<var_dim;i++) var_label[i] = "~";
  21.     /* secondary phase space variable label (DIM=var_dim) */
  22.     for(i=0;i<var_dim;i++) var_polar_label[i] = "~";
  23.     /* parameter variable label (DIM=param_dim)*/
  24.     for(i=0;i<param_dim;i++) param_label[i] = "~";
  25.     /* function variable label (DIM=func_dim)*/
  26.     for(i=0;i<func_dim;i++) func_label[i] = "~";
  27.     /* starting parameter values (DIM=param_dim)*/
  28.     for(i=0;i<param_dim;i++) param[i] = 0;
  29.     /* starting primary phase space variable values (DIM=var_dim)*/
  30.     for(i=0;i<var_dim;i++) var_i[i] = 0;
  31.     /* starting seconsary phase space variable values (DIM=var_dim)*/
  32.     for(i=0;i<var_dim;i++) var_polar_i[i] = 0;
  33.     /* starting bounds of parameter window box */
  34.     for(i=0;i<param_dim;i++) {
  35.         param_min[i] = -5;
  36.         param_max[i] = 5;
  37.     }
  38.     /* starting bounds of primary phase space window box */
  39.     for(i=0;i<var_dim;i++) {
  40.         var_min[i] = -5;
  41.         var_max[i] = 5;
  42.     }
  43.     /* starting bounds of secondary phase space window box */
  44.     for(i=0;i<var_dim;i++) {
  45.         var_polar_min[i] = -5;
  46.         var_polar_max[i] = 5;
  47.     }
  48.  
  49.     /* period */
  50.         if(mapping_on){
  51.         ir = 1;
  52.     }
  53.     else {
  54.         ir = 0;
  55.     }
  56.  
  57.     /* variable projection plane */
  58.     var_cur_index[0] = 0;
  59.     var_cur_index[1] = 1; 
  60.  
  61.     /* parameter projection plane */
  62.     param_cur_index[0] = 0; 
  63.     param_cur_index[1] = 1; 
  64.  
  65.     /* offset for a periodic orbit */
  66.     for(i=0;i<var_dim;i++)
  67.         xoffset[i]=0;
  68.  
  69.     /* rescaling factor for QC integrators */
  70.     for(i=0;i<var_dim;i++)
  71.         int_yscal[i]=1.;
  72.  
  73.     /* auxiliary window parameters */
  74.     for(i=0;i<aux_max;i++){
  75.         aux_on[i] = 0;
  76.         aux_win_mode[i] = 0;
  77.         aux_rescale[i] = 0;
  78.         aux_cur_index[i][0] = 0;
  79.         aux_cur_index[i][1] = 1;
  80.     }
  81. }
  82.